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Abstract 

We compute explicit upper bounds on the distance between the law of a multivari- 
ate Gaussian distribution and the joint law of wavelets/needlets coefficients based on 
a homogeneous spherical Poisson field. In particular, we develop some results from 
Peccati and Zheng (2011), based on Malliavin calculus and Stein's methods, to assess 
the rate of convergence to Gaussianity for a triangular array of needlet coefficients 
with growing dimensions. Our results are motivated by astrophysical and cosmolog- 
ical applications, in particular related to the search for point sources in Cosmic Rays 
data. 
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1 Introduction 



The aim of this paper is to establish multidimensional normal approximation results for 
vectors of random variables having the form of wavelet coefficients integrated with respect 
to a Poisson measure on the unit sphere. The specificity of our analysis is that we require 
the dimension of such vectors to grow to infinity. Our techniques are based on recently 
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obtained bounds for the normal approximation of functionals of general Poisson measures 
(see [HJ S3]), as well as on the use of the localization properties of wavelets systems on 
the sphere (see [35], as well as the recent monograph [30J). A large part of the paper is 
devoted to the explicit determination of the above quoted bounds in terms of dimension. 

1.1 Motivation and overview 

A classical problem in asymptotic statistics is the assessment of the speed of convergence 
to Gaussianity (that is, the computation of explicit Berry-Esseen bounds) for parametric 
and nonparametric estimation procedures - for recent references connected to the main 
topic of the present paper, see for instance [T7J [29J, [55]. In this area, an important novel 
development is given by the derivation of effective Berry-Esseen bounds by means of the 
combination of two probabilistic techniques, namely the Malliavin calculus of variations 
and the Stein's method for probabilistic approximations. The monograph [8] is the standard 
modern reference for Stein's method, whereas [38] provides an exhaustive discussion of the 
use of Malliavin calculus for proving normal approximation results on a Gaussian space. 
The fact that one can use Malliavin calculus to deduce normal approximation bounds (in 
total variation) for functionals of Gaussian fields was first exploited in [37J - where one 
can find several quantitative versions of the "fourth moment theorem" for chaotic random 
variables proved in [39]. Lower bounds can also be computed, entailing that the rates of 
convergence provided by these techniques are sharp in many instances - see again [38J. 

In a recent series of contributions, the interaction between Stein's method and Malliavin 
calculus has been further exploited for dealing with the normal approximation of functionals 
of a general Poisson random measure. The most general abstract results appear in [32J (for 
one-dimensional normal approximations) and [13] (for normal approximations in arbitrary 
dimensions). These findings have recently found a wide range of applications in the field 
of stochastic geometry - see [251 ESI E31 [23 SB] for a sample of geometric applications, as 
well as the webpage 

http : / /www . iecn . u-nancy . f r/ ~nourdin/steinmalliavin . htm 

for a constantly updated resource on the subject. 

The purpose of this paper is to apply and extend the main findings of [HI 03] in order to 
study the multidimensional normal approximation of the elements of the first Wiener chaos 
of a given Poisson measure. Our main goal is to deduce bounds that are well-adapted to 
deal with applications where the dimension of a given statistic increases with the number of 
observations. This is a framework which arises naturally in many relevant fields of modern 
statistical analysis; in particular, our principal motivation originates from the implemen- 
tation of wavelet systems on the sphere. In these circumstances, when more and more data 
become available, a higher number of wavelet coefficients is evaluated, as it is customarily 
the case when considering, for instance, thresholding nonparametric estimators. We shall 
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hence be concerned with sequences of Poisson fields, whose intensity grows monotonically. 
We then exploit the wavelets localization properties to establish bounds that grow linearly 
with the number of functionals considered; we are then able to provide explicit recipes, for 
instance, for the number of joint testing procedures that can be simultaneously entertained 
ensuring that the Gaussian approximation may still be shown to hold, in a suitable sense. 



1.2 Main contributions 

Consider a sequence {Xj : i > 1} with values in the unit sphere S 2 , and define {4>jk} 
to be the collection of the spherical needlets associated with a certain constant B > 1, 
see Section I3TT1 below for more details and discussion. Write also a^ k = E[ipj k (Xi) 2 } and 
bjk = E[i[)jk(Xi)], and consider an independent (possibly inhomogeneous) Poisson process 
{N t : t > 0} on the real line such that E[N t ] = R(t) — > oo, as t — > oo. Formally, our 
principal aim is to establish conditions on the sequences {j{n) : n > 1}, {R{n) : n > 1} 
and {d(n) : n > 1} ensuring that the distribution of the centered <i(n)-dimensional vector 

Y n = (Y n ,li Yn,d{n)) (1-1) 

^(n)b {Xj) _ R{n)b j{n)kl ^ ^j(n)k d(n) {Xi) R{n)b j{n)kd{n) 



\/R{n) \ ~[ ^i(n)fei Pj{n)ki ' ^ cr i(n)fc 1 Vj(n)k, 



d{n) 



is asymptotically close, in the sense of some smooth distance denoted d<i (see Definition 
12. 6p . to the law of a <i(n)-dimensional Gaussian vector, say Z n , with centered and indepen- 
dent components having unit variance. The use of a smooth distance allows one to deduce 
minimal conditions for this kind of asymptotic Gaussianity. The crucial point is that we 
allow the dimension d(n) to grow to infinity, so that our results require to explicitly assess 
the dependence of each bound on the dimension. We shall perform our tasks through 
the following main steps: (i) Proposition 14.11 deals with one-dimensional normal approxi- 
mations, (ii) Proposition 15.41 deals with normal approximations in a fixed dimension, and 
finally (iii) in Theorem 15.51 we deduce a bound that is well-adapted to the case d{n) — > oo. 
More precisely, Theorem 15.51 contains an upper bound linear in d(n), that is, an estimate 
of the type 

d 2 (Y n , Z n ) < C(n) x d(n) . (1.2) 



It will be shown in Corollary 15. 6[ that the sequence C(n) can be chosen to be 



O ll/^R{n)B- 2 ^) 



as discussed below in Remark I4.3[ R(n) x B 2j(n ) can be viewed as a measure of the 
"effective sample size" for the components of Y n . 
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1.3 About de-Poissonization 



Our results can be used in order to deduce the asymptotic normality of de-Poissonized 
linear statistics with growing dimension. To illustrate this point, assume that the random 
variables Xi are uniformly distributed on the sphere. Then, it is well known that bjk = 0, 
whenever j > 1. In this framework, when j[n) > 1 for every n, R(n) = n and d{n)/n l / A — y 

0, the conditions implying that Y n is asymptotically close to Gaussian, automatically ensure 
that the law of the de-Poissonized vector 

V™ o- j{n)kl a j{n)kd(n) J 

is also asymptotically close to Gaussian. The reason for this phenomenon is nested in the 
statement of the forthcoming (elementary) Lemma 11.11 

Lemma 1.1 Assume that R(n) = n, that the Xi 's are uniformly distributed on the sphere, 
and that j(n) > 1 for every n. Then, there exists a universal constant M such that, for 
every n and every Lipschitz function ip : IR d(n ) — > WL, the following estimate holds: 

din) 

E[<p{Y$] - E[ v {Y n )} < M|MU ip -^i 
Proof. Fix I = 1, d(n), and write (3i{x) = ^ J(n)fc ' — , in such a way that E[(3i(Xi) 2 ] = 

j(n)ki 

1. One has that 

E[(Y^ :l -Y n j) 2 ]= 2(1 -a n ), 

where 

1 e~ n n m e~ n n n 

a n — — y ; — (n Am) = 1 ; — . 

n ^— ' 

m=0 

This gives the estimate 



n c — ' ml ni 

m=0 



ElKl - Yn,l\] < y/ElKl - Yn,?} < V 2 ^T' 

so that the conclusion follows from an application of Stirling's formula and of the Lipschitz 
property of if. ■ 

Remark 1.2 (i) Lemma implies that one can obtain an inequality similar to M.2fy 
for Y' n , that is: 

d 2 (Y;,Z n )< fc(n) + ^\ xd(n). 
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(ii) With some extra work, one can obtain estimates similar to those in Lemma \l.l\ also 
when the constants 6j( n )fc i are possibly different from zero. This point, that requires 
some lenghty technical considerations, falls slightly outside the scope of this paper and 
will be pursued in full generality elsewhere. 

(Hi) In Jflj, Bentkus proved the following (yet unsurpassed) bound. Assume that {Xi : i > 
1} is a collection of i.i.d. d-dimensional vectors, such that X\ is centered and with 
covariance equal to the identity matrix. Set S n = nr 1 ' 2 {X\ + ■ ■ ■ X n ), n > 1 and let 
Z be a d-dimensional centered Gaussian vector with i.i.d. components having unit 
variance. Then, for every convex set C C M. d 



E[lc(S n )} - E[1 C (Z)\ 



< 
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where (3 = E[\\Xi\\^ d ]. It is unclear whether one can effectively use this bound in 
order to investigate the asymptotic Gaussianity of sequences of random vectors of 
the type U. ip -z TOj) . in particular because, for a fixed n, the components of Y n , 
have in general a non trivial correlation. Note also that a simple application of 
Jensen inequality shows that (3d 1 ' /4 n -1 / 2 > d 1 l i n~ x l 2 . However, a direct comparison 
of Bentkus' estimates with our "linear" rate in d (see M.2) . as well as Theorem \5.5\ 
below) is unfeasible, due to the differences with our setting, namely concerning the 
choice of distance, the structure of the considered covariance matrices, the Poissonized 
environment, and the role of B^ n ^ discussed in Remark\4- 3 . 



(iv) A careful inspection of the proofs of our main results reveals that the findings of this 
paper have a much more general validity, and in particular can be extended to kernel 
estimators on compact spaces satisfying mild concentration and equispacing properties 
(see also I2ty 121)/). In this paper, however, we decided to stick to the presentation 
on the sphere for definiteness, and to make the connection with applications clearer. 
Some more general frameworks are discussed briefly at the end of Section^ 

(v) For notational simplicity, throughout this paper we will stick to the case where all the 
components in our vector statistics are evaluated at the same scale j(n) (see below for 
more precise definitions and detailed discussion). The relaxation of this assumption 
to cover multiple scales (j±(n), ...j d {n)) does not require any new ideas and is not 
considered here for brevity's sake. 



1.4 Plan 

The plan of the paper is as follows: in Section [2] we provide some background material on 
Stein-Malliavin bounds in the case of Poisson random fields, and we describe a suitable 
setting for the current paper, entailing sequences of fields with monotonically increasing 
governing measures. We provide also some new results, ensuring that the Central Limit 
Theorems we are going to establish are stable, in the classical sense. In Section [3] we recall 
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some background material on the construction of tight wavelet systems on the sphere (see 
[5^135] for the original references, as well as [30, Chapter 10]) and we explain how to express 
the corresponding wavelet coefficients in terms of stochastic integrals with respect to a 
Poisson random measure. We also illustrate shortly some possible statistical applications. 
In Section H] we provide our bounds in the one- dimensional case; these are simple results 
which could have been established by many alternative techniques, but still they provide 
some interesting insights into the "effective area of influence" of a single component of the 
wavelet system. The core of the paper is in Section where the bound is provided in the 
multidimensional case, allowing in particular for the number of coefficients to be evaluated 
to grow with the number of observations. This result requires a careful evaluation of the 
upper bound, which is made possible by the localization properties in real space of the 
wavelet construction. 

2 Poisson Random Measures and Stein-Malliavin 
Bounds 

In order to study the asymptotic behaviour of linear functionals of Poisson measures on 
the sphere § 2 , we start by recalling the definition of a Poisson random measure - for more 
details, see for instance jH] H7J [50]. We work on a probability space (ft, J 7 , P). 

Definition 2.1 Let [Q,A, /x) be a a-finite measure space, and assume that fi has no atoms 
(that is, fi({x}) = 0, for every x G Q). A collection of random variables {N (A) : A G A} , 
taking values in Z + U {+00} , is called a Poisson random measure (PRM) on G with 
intensity measure (or control measure) /i if the following two properties hold: 

1. For every A G A, N (A) has Poisson distribution with mean fi (A); 

2. If Ai, . . . A n G A are pairwise disjoint, then N (Ai) , . . . , N (A n ) are independent. 

Remark 2.2 (i) In Definition \2.1\ a Poisson random variable with parameter A = 00 
is implicitly set to be equal to 00. 

(ii) Points 1 and 2 in Definition \2.1\ imply that, for every u G Vt, the mapping A 1— > 
N(A,u) is a measure on 0. Moreover, since fi is non atomic, one has that 

P[N({x}) = or 1, Vx G 9] = 1. (2.4) 



Assumption 2.3 Our framework for the rest of the paper will be the following special case 
of Definition \2.1l 

(a) We take 6 = M + x S 2 , with A = £>(©), the class of Borel subsets of 0. 
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(b) The symbol N indicates a Poisson random measure on B, with homogeneous intensity 
given by p. = p x v, where p is some measure on ¥L + and v is a 'probability on § 2 of 
the form v(dx) = f(x)dx, where f is a density on the sphere. We shall assume that 
P({0}) = an d ^at the mapping p i— > p([0,t\) is strictly increasing and diverging to 
infinity as t — > oo. We also adopt the notation 

Rt:=p([0,t]), t>0, (2.5) 

that is, t R t is the distribution function of p. 

Remark 2.4 

(%) For a fixed t > 0, the mapping 

A t-> N t (A) := N([0, t] x A) (2.6) 

defines a Poisson random measure on § 2 , with non-atomic intensity 

Pt(dx) = R t ■ v(dx) = R t ■ f(x)dx. (2.7) 

Throughout this paper, we shall assume f(x) to be bounded and bounded away from 
zero, e.g. 

Ci < /O) < C2 , some Ci, C2 > , for all x G § 2 . (2.8) 

(ii) Let {Xi = i > 1} be a sequence of i.i.d. random variables with values in S 2 and 
common distribution equal to v. Then, for a fixed t > 0, the random measure 
A 1 — Y N t (A) = N([0,t] x A) has the same distribution as A ^ ^2 i=1 6xi(A), were 
8 X indicates a Dirac mass at x, and N is an independent Poisson random variable 
with parameter R t . This holds because: (a) since v is a probability measure, the 
support of the random measure N t (written supp(N t )) is almost surely a finite set, 
and (b) conditionally on the event {N t (S 2 ) = n} (which is the same as the event 
{ \ supp(N t )\ — n} - recall l[2.4\ )), the points in the support of N t are distributed as n 
i.i.d. random variables with common distribution v. 

(Hi) By definition, for every t\ < t 2 one has that a random variable of the type N t2 (A) — 
N tl (A), A C § 2 , is independent of the random measure N tl , as defined in 112.6}) . 

(iv) To simplify the discussion, one can assume that p(ds) = R ■ £(ds), where i is the 
Lebesgue measure and R > 0, in such a way that R t = R ■ t. 

We will now introduce two distances between laws of random variables taking values 
in ¥L d . Both distances define topologies, over the class of probability distributions on M. d , 
that are strictly stronger than convergence in law. One should observe that, in this paper, 



7 



the first one (Wasserstein distance) will be only used for random elements with values in 
E. Given a function g E C 1 (IR d ), we write 1 1 1 1 iip = sup || Vg(x)\\ R d. If g E C 2 (R d ), we set 

M 2 (g) = sup ||Hess^(x)||op, 

x€R d 

where || • || op indicates the operator norm. 

Definition 2.5 The Wasserstein distance dw, between the laws of two random vectors 
X,Y with values inM. d (d>l) and such that E ||X|| Rd ,E \\Y\\ Rd < oo, is given by: 

d w (X, Y) = sup \E [g (X)] — E[g (Y)} \ , 
g-\\9\\Li P <t 

Definition 2.6 The distance d 2 between the laws of two random vectors X, Y with values 
in M. d (d > 1 ), such that E \\X\\ Rd , E ||^|| R d < oo ; is given by: 

d 2 (X,Y)=sup\E[g(X)]-E[g(Y)}\, 
gen 

where H denotes the collection of all functions g E C 2 (R d ) such that \\g\\up ^ 1 an d 
M 2 (g) < 1. 

We now present, in a form adapted to our goals, two upper bounds involving random 
variables living in the so-called first Wiener chaos of N. The first bound was proved in 
|42j . and concerns normal approximations in dimension 1 with respect to the Wasserstein 
distance. The second bound appears in [13], and provides estimates for multidimensional 
normal approximations with respect to the distance d 2 . Both bounds are obtained by 
means of a combination of the Malliavin calculus of variations and the Stein's method for 
probabilistic approximations. 

Remark 2.7 (i) Let f E L 2 (Q,n) flL 1 (@,/i). In what follows, we shall use the symbols 
N(f) and N(f), respectively, to denote the Wiener-Ito integrals of f with respect to 
N and with respect to the compensated Poisson measure 

N(A) = N(A) - n(A), AeB(Q), (2.9) 

where one uses the convention N(A) — fi(A) = oo whenever fi(A) = oo (recall that 
fi is a -finite). Note that, for N(f) to be well-defined, one needs that f E L 1 (0,/^) ; 
whereas for N(f) to be well-defined one needs that f E L 2 {Q,n). We will also make 
use of the following isometric property: for every f,g EL 2 (Q, fi) , 

E[N{f)N{g)} = [ f(x)g(x)fx(dx). (2.10) 

The reader is referred e.g. to fjlj Chapter 5] for an introduction to Wiener-Ito 
integrals. 
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(ii) For most of this paper, we shall consider Wiener-Ito integrals of functions f having 
the form f = [0, t] x h, where t > and h G L 2 (S> 2 , v) fl -^(S 2 , v). For a function f 
of this type one simply writes 

N(f) = N([0, t] x h) := N t (h), and N(f) = N([0, t] x h) := N t (h). (2.11) 



Observe that this notation is consistent with the one introduced in Ii2.6\) . Indeed, it 
is easily seen that N t (h) (resp. N t (h)) coincide with the Wiener-Ito integral of h 
with respect to N t (resp. with respect to the compensated measure N t — N t — m — 
N t -Rf v). 



(Hi) In view of Remark\2.4\-(H), one also has that, for h G L (§ , v) fl L (S> , v), 



N t (h) 



Y] h{x), andN t (h) = V h{x)- I 

, at \ ^ , at ^ 



h(x)n t (dx), (2.12) 



xGsupp(Nt) 



xGsupp(Nt) 

with fi t defined as in \2. 7| ). 

Theorem 2.8 Let the notation and assumptions of this section prevail. 

1. Let h G L 2 (S 2 , v) := L?{y), let Z ~ ^(0, 1) andfixt>0. Then, the following bound 
holds (remember the definition \2. J\j): 



d w (N t (h),Z) < 



1 



£ 2 (§Vt) 



+ / \h(z)\^ t (dz). 

§ 2 



(2.13) 



As a consequence, if {h t } C L 2 [y) fl L 3 (u) is a collection of kernels verifying, as 
t — > oo, 



II^IU 2 (s 2 ,Mt) ~^ 1 and ll^ilU 3 (§ 2 ,M t ) -> °> 

one /ias t/ie CLT 



N(ht) ^ 



(2.14) 



(2.15) 



and the inequality \2.13) provides an explicit upper bound in the Wasserstein distance. 
2. For a fixed integer d > 1, letY~ jV d (0, C), with C positive definite and let 

F t = (F tA , . . . , F t , d ) = (N t (ht,x) ,...N t (Kd) 

be a collection of d- dimensional random vectors such that h tia G L?(y). If we call T t 
the covariance matrix of F t , that is, 

T t (a,b) = E N t (ht, a ) N t (ht,b) = (h t , a , Kb) l 2 (s 2 ,^) > a > b = 1 ' -> d ' 
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then: 



d 2 (F t ,Y) < \\C-\JC\\l\\C-T t \\ H ^ 



^7T 



where 



(2.16) 

C_1 |l™ II^IL E / \ht,i(x)\\h t:j (x)\\h tik (x)\^ t (dx) , 
l|C|ltll Cr - r *Hw. (2- 17 ) 

d „ 

c -1 III llciU E / I^WIVW, 

and || • ||h.s. stand, respectively, for the operator and Hilbert- Schmidt 



< He- 1 

4 



op 11 M op 



norms. In particular, if T t (a, 6) — > C (a, 6) and L 2 \ht, a (x)\ fit (dx) — > as t - 
oo, for a,b = 1, . . . d, then d 2 (F t , Y) — > and F t converges in distribution to Y. 



Remark 2.9 The estimate 112.16]) will be used to deduce one of the main multidimensional 
bounds in the present paper. It is a direct consequence of Theorem 3.3 in f4^i , where the 
following relation is proved: for every vector (Fi,...,F d ) of sufficiently regular centered 
functionals of N t , 



d 2 (F,X)<\\C^\\jC\\l 



a 



Oh) 



+ ^ ||c _ 1||3 / 2||c| 



I op 



op 



fj,t(dz)K 



Ei^i) (El^ 1 ^ 

3=1 



i=l 



where 



and 



D z F(uj) = F z (u) - F(u) , a.e. - fi(dz)P(dco) , 



F Z (N) = F Z {N + 6 Z ), 



that is, the random variable F z is obtained by adding to the argument of F (which is a 
function of the point measure N ), a Dirac mass at z, and L~ x is the so-called pseudo-inverse 
of the Ornstein- Uhlenbeck operator. The estimate 112.16}) is then obtained by observing that, 
when Fi = F ti = N t (h tti ), then D z Fi = —D Z L~ 1 F = h tti (z), in such a way that 



\ 



a 



\ C ~ K t\\ H .S. ' 



'■J 
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and 



Ei^i [12\ D * L ~ lF i 

j=i j \j=i 

d 

£ 

i,j,k=l 



I h t ,i (x) | | h tJ (x) | | h t ,k (x) | fit (dx) . 



The next statement deals with the interesting fact that the convergence in law implied 
by Theorem 12.81 is indeed stable, as defined e.g. in the classic reference [19, Chapter 4]. 

Proposition 2.10 The central limit theorem described at the end of Point 2 of Theorem 
\2.8\ (and a fortiori the CLT at Point 1 of the same theorem) is stable with respect to o~(N) 
(the a- field generated by N) in the following sense: for every random variable X that is 
a (N) -measurable, one has that 

(X,F t )^(X,Y), 

where Y ~ ^V d (0,C) is independent of N. 

Proof. We just deal with the case d — 1, the extension to a general d following 
from elementary considerations. A density argument shows that it is enough to prove the 
following claim: if N(h n ) (h n G L 2 (fi), n > 1) is a sequence of random variables verifying 
E[N(h n ) 2 ] = \\h n \\ 2 L 2^ — > 1 and J e \h n \ 3 dfi — > 0, then for every fixed / G L 2 (fi), the pair 

(N(f), N(h n )) converges in distribution, as n — > oo, to (N(f),Z), where Z ~ r yV(0,l) is 
independent of N. To see this, we start with the explicit formula (see e.g. [4"Tj formula 
(5.3.31)]): for every A, 7 G R 

M\l)-= E[exp(t\N(f)+lN(h n ))} 



exp 



,iXf(x)+i'yh„(x) 



l — i(\f(x) + , yh n (x))] n(dx 



e 



Our aim is to prove that, under the stated assumptions, 



lim log(V>»(A,7)) 



e 



r 
2 



Standard computations show that 

r 2 
log(^ n (A, 7 )) - { J[e lXnx) - l-i\f(x)] fi(dx) - X-} 



< 



7 2 7 2 f l7| 3 f 

y - y J K{x) n(dx) + I7AI \ (h n ,f) L 2 {fl) \ + —J \h n (x)\ fi(dx) . 



Since J e \h n (x)\ 3 ^(dx) —> and the mapping n \-> H^nll^^) is bounded, one has that 
{h n ,f)L 2 (fi) 0, and the conclusion follows by using the fact that H/inJIla^) - > 1 by 
assumption. ■ 
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3 Needlet coefficients 



3.1 Background: the needlet construction 

We now provide an overview of the construction of the set of needlets on the unit sphere. 
The reader is referred to [301 Chapter 10] for an introduction to this topic. Relevant 
references on this subject are: the seminal papers [351 |36J, where needlets have been 
first defined; [131 EH CE21 EES], among others, for generalizations to homogeneous spaces of 
compact groups and spin fiber bundles; [31 HI [2B1 132] for the analysis of needlets on spherical 
Gaussian fields, and [311 SHI El EH] for some (among many) applications to cosmological and 
astrophysical issues; see also [3HH9] for other approaches to spherical wavelets construction. 

(Spherical harmonics) In Fourier analysis, the set of spherical harmonics 

{Y lm :l> 0,m = -/,..., 1} 

provides an orthonormal basis for the space of square-integrable functions on the unit 
sphere L 2 (S 2 ,dx) := L 2 (§ 2 ), where dx stands for the Lebesgue measure on § 2 (see for 
instance [21 1231 EH1 EI])- Spherical harmonics are defined as the eigenfunctions of the 
spherical Laplacian A S 2 corresponding to eigenvalues — 1(1 + 1), e.g. A S 2Yi m = —1(1 + 
l)Yi m , see again [301 EH EI] for analytic expressions and more details and properties. For 
every / > 0, we define /Q as the linear space given by the restriction to the sphere of the 
polynomials with degree at most /. Plainly, one has that 

i 

^ = (B span { Ykm : m = ~ k > ^ ' 

fc=0 

where the direct sum is in the sense of L 2 (§> 2 ). 

(Cubature points) It is well-known that for every integer I = 1, 2, ... there exists a finite set 
of cubature points Qi C § 2 , as well as a collection of weights { A^}, indexed by the elements 
of Qi, such that 

V/ e K h [ f(x)dx = V X v f( v ). 

Now fix B > 1, and write [x] to indicate the integer part of a given real x. In what follows, 
we shall denote by Xj = {C,jk} and {Xjk}, respectively, the set Qpsj+i] and the associated 
class of weights. We also write Kj = card{Xj}. As proved in (35J, [36] , cubature points and 
weights can be chosen to satisfy 

X jk w B- 21 , Kj ■« B 21 , (3.18) 

where by a b, we mean that there exists ci, c% > such that c±a <b< c%a (see also e.g. 
[3 m HE] and [IS Chapter 10]). 
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(Spherical needlets) Fix B > 1 as before, as well as a real- valued mapping b on (0, oo). 
We assume that b verifies the following properties: (i) the function b (•) has compact 
support in [B^ 1 , B] (in such a way that the mapping I \— > b (-^j) has compact support in 
/ G [-B- 7-1 , (ii) for every £ > 1, Y^jLob 2 (£B~i) = 1 (partition of unit property), and 

(hi) b(-) G C°° (0, oo). The collection of spherical needlets {ipjk}, associated with B and 
&(•), are then defined as a weighted convolution of the projection operator Li((x,y)) = 
Y! m =-i Y im (%) Yirn (y) , that is 

ipjk(x) := \^^2 b y^j L i(( x ^jk)) , (3.19) 



(Localization) The properties of b entail the following quasi-exponential localization prop- 
erty (see [35J or [30, Section 13.3]): for any r = 1,2,... there exists k t > such that for 
any x G § 2 , 

£ (i + gJ^kfa))r (3 - 20) 

where d(x,y) := arccos ((x, y)) is the spherical distance. From localization, the following 
bound can be established on the L p (S 2 ) norms: for all 1 < p < +oo, there exist two 
positive constants q p and q' p such that 

q p B>M ] < Ujk\\ Lp{s2) < q'M 1 -!). (3.21) 

(Needlets as frames) Finally, the fact that b is a partition of unit, allows on to deduce the 
following reconstruction formula (see again [35]): for / G L 2 (§ 2 ): 

f( x ) = ^2l3jk*Pjk(x) , 

where the convergence of the series is in L 2 (§ 2 ), and 

Pjk (f, ^jk) L2(S 2) = / / (x) ip jk (x) dx , (3.22) 

Js 2 

represents the so-called needlet coefficient of index j, k. 

3.2 Two motivations: density estimates and point sources 

The principal aim of this paper is to establish multidimensional asymptotic results for some 
possibly randomized version of random variables of the type 

1 n 

?J* = ®? = -$>J-*M' j = k = l,...,Kj, (3.23) 



n 

i=l 
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where the function ipj k is defined according to (13.191) . and {Xj : % > 1} is some adequate 
sequence of i.i.d. random variables. We may also study the asymptotic behaviour, as 

t — > oo, of mult i- dimensional object of the type k = 1, 2, Kj(t) \, where t i— > Kj(t) 

is a non- decreasing mapping possibly diverging to infinity, and j may change with t. In 
other words, as happens in realistic experimental circumstances, we may decide to focus on 
a growing number of coefficients as the number of (expected) events increase. Two strong 
motivations for this analysis, both coming from statistical applications, are detailed below. 

(Density estimates) Consider a density function / on the sphere S 2 , that is: / is a 
mapping from § 2 into R + , verifying J g2 f(x)dx = 1, where dx indicates the Lebesgue 
measure on § 2 . Let {Xi : % = 1, ...,n} be a collection of i.i.d. observations with values in 
S> 2 with common distribution given by f(x)dx. A classical statistical problem, considered 
for instance by E2J, [25j , concerns the estimation of / by wavelets/needlets thresholding 
techniques. To this aim, keeping in mind the notation (I3.23p . one uses ([ID]. [16]) the 
following estimator of /: 

fa) = E (x) , Pf k := W { |^|> ctn }, 

jk 

where t n = A/log n/n and c is a constant to be determined. Finite-sample approximations 
on the distributions of (3jk can then be instrumental for the exact determination of the 
thresholding value ct n , see e.g. [TO] IT6]. 

(Searching for point sources) The joint distribution of the coefficients {Pjk} (as defined 
in (I3.23P ) is required in statistical procedures devised for the research of so-called point 
sources, again for instance in an astrophysical context (see for instance [53]). The physical 
issue can be formalized as follows: 

- Under the null hypothesis, we are observing a background of cosmic rays governed 
by a Poisson measure on the sphere § 2 , with the form of the measure N t (-) defined in 
( 12.61) for some t > 0. In particular, N t is built from a measure N verifying Assumption 
(12.31) . and the intensity of Ht(dx) = E[N t (dx)] is given by the absolutely continuous 
measure R t ■ f(x)dx, where R t > and / is a density on the sphere. This situation 
corresponds, for instance, to the presence of a diffuse background of cosmological 
emissions. 

- Under the alternative hypothesis, the background of cosmic rays is generated by a 
Poisson point measure of the type: 




where {£i,...,£p} C S 2 , each mapping t h-> Nj: p ' is an independent Poisson process 
over [0, oo) with intensity X p , and 

yjc,{z)dz = l\ e A] . 
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In this case, one has that N£ is a Poisson measure with atomic intensity 




In this context, the informal expression "searching for point sources" can then be trans- 
lated into "testing for P = 0" or "jointly testing for X p > at p = 1, ...,P". The number 
P and the locations {£i, ...£p} can be in general known or unknown. We refer to [TSl E2] 
for astrophysical applications of these ideas. 

Remark 3.1 In order to directly apply the findings of f4^ \4^j , in what follows we shall 
focus on a randomized version of A3. 23\) . where n is replaced by an independent Poisson 
number whose parameter diverges to infinity. Also, we will prefer a deterministic normal- 
ization over a random one. As formally shown in the discussion to follow, the resulting 
randomized coefficients can be neatly put into the framework of Section 

3.3 Needlet coefficients as Wiener-Ito integrals 

Let N be a Poisson measure on R + x § 2 satisfying the requirements of Assumption 12.31 
(in particular, the intensity of N has the form p x i/, where u(dx) = f(x)dx, for some 
probability density / on the sphere, and one writes R t = p([0,t]), t > 0). For every 
t > 0, let the Poisson measure iV* on S 2 be defined as in (12. 6p . For every j > 1 and every 
k = l,...,Nj, consider the function ipjj. defined in (I3.19p . and observe that ipjk is trivially 
an element of L 3 (§ 2 , v) n L 2 (§ 2 , v) n L 1 ^ 2 , v). We write 

Observe that, if f(x) = j- (that is, the uniform density on the sphere), then bjk = for 
every j > 1. On the other hand, under (12. 8p . 

CiU j k(.)\\l 2 <CT%<C2U j k(.)\\l 2 ■ (3.24) 

Note that (see I3.2ip the L 2 -norm of {ipjk} is uniformly bounded above and below, and 
therefore the same is true for {cr 2 fc j (indeed, there exists n > 0, independent of j and k, 
such that < n < \\^jk\\ 2 L 2t§2\ < !)• For every t > and every j, k, we introduce the kernel 
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In view of Remark l2.4j -(ii). the random variable p^f* can always be represented in the form 



(Rt) 



t<Jjk 



where {JQ : i > 1} is a sequence of i.i.d. random variables with common distribution u, 
and independent of the Poisson random variable N t (E> 2 ). Moreover, the following relations 
are immediately checked: 



Remark 3.2 Using the notation $3.23\) . we have that 
Pjk ~ 



(3.27) 



4 Bounds in dimension one 

We are now going to apply the content of Theorem 12.81 (1) to the random variables (3^^ 
introduced in the previous section. In the next statement, we write Z ~ A/"(0, 1) to indicate 
a centered Gaussian random variable with unit variance. Recall that £ 2 : = sup^es 2 |/ (x)\, 
p > 1, and that the constants q p , q' p have been defined in (I3.2ip . 



Proposition 4.1 For every j, k and every t > 0, one has that 



R 



It follows that for any sequence (j(n), k(n),t(n)) } Pjf n )k(n) conver 9 es ^ n distribution to Z, 
as n oo, provided B 2 ^'^ = o(R t ( n )). The convergence is a (N) -stable, in the sense of 
Proposition \2.1(A 

Proof. Using fl3T2oD - fl3T2oP together with fl2TT7j) and (J23D, 



d w {^\z) < / 



hf ] (x) 



< 



VR]" 3 
VRto-% 



fit (dx) 
\4>jk (x)\ 3 f(x) dx < 



l R~t (J % 



J jk\\L 3 (S 2 ) 



where in the last inequality we use the property (I3.2ip with p = 3 to have: 



< (q' 3 fBVM) = {q',fBK 
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The last part of the statement follows from the fact that the topology induced by the 
Wasserstein distance (on the class of probability distributions on the real line) is strictly 
stronger than the topology of convergence in law. ■ 



Remark 4.2 For f(x) = {47r} we have 

°% = ^ J §2 rfk ( x ) dx = WjkWvQP) > 
and more generally, under K2. 8\) . 

(Rt) „\ . & (^) 3 C: 



d W (p\?>,Z)<-= ^ - := 7 (j,M) . (4.28) 

Remark 4.3 The previous result can be given the following heuristic interpretation. The 
factor B^i can be viewed as the "effective scale" of the wavelet, i.e. it is the radius of 
the region centred at £jk where the value of the wavelet function is not negligible. Because 
needlets are isotropic, the "effective area" is of order B~ 2 K For governing measures with 
density which is bounded and bounded away from zero, the expected number of observations 
on a spherical cap of radius i?"- 7 around £jk is hence given by 



E [card {X t : d(X h & fc ) < B~i}] ~ R t I f{x)dx , 

J d(x£ jk )<B-j 

(iB- 2j R t <R t f f(x)dx < ( 2 B- 2j R t , 



1 d(x,£ jk )<B-i 

using U\ equation (8)]. Because the Central Limit Theorem can hold only when the effective 
number of observations grows to infinity, the condition B~ 2 ^R t — > oo is quite expected. In 
the thresholding literature, coefficients are usually considered up to the frequency Jr such 
that B 2Jr ~ R t /\ogR t , see for instance fl6] / and f^/; under these circumstances, we have 



1 



d 2 [f^lz] = O ( -^== )^0forR t ^ +oo 



Therefore /3j^j does converge in law to Z. 



5 Multidimensional bounds 

We are now going to apply Part 2 of Theorem (12. 8 p to the computation of multidimensional 
Berry-Esseen bounds involving vectors of needlet coefficients of the type (I3.26p . After 
having proved some technical estimates in Section 15.14 we will consider two bounds. One 
is proved in Section l5~2l by means of (I2.17p . and it is well adapted to the case where the 
number of needlet coefficients, say d, is fixed. In Section [5~3| we shall focus on (I2.16p . and 
deduce a bound which is adapted to the case where the number d is possibly growing to 
infinity. 
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5.1 A technical result 



The following estimate, allowing one to bound the covariance between any two needlet 
coefficients, will be used throughout this section. We let the notation and assumptions of 
the previous section prevail. 

Lemma 5.1 For any j > 1 and k\ ^ k 2 < Kj = card{Xj} and every t > 0, there exists 
a constant C T > 0, solely depending on r, and such that 



\T Rt (h,k 2 )\:=\E^^\< 



C T ( 2 



o~jk 1 o~jk 2 (1 + Bid (£,jk 1 ,£,jk 2 )) T 



Proof. We focus on r > 2; note that the inequality for any fixed value of r immediately 
implies the result for all r' < r. For k\ ^ k 2 we have: 

\Y Rt (k u k 2 )\ = 



< 



Rt& j kl 0~ jk 2 

Rt 

RtO'jk 1 0~jk 2 
(2 



■0i*i ( x ) ipjki ( x ) V>t{dx 
ip jkl (x) ip jk2 (x) f(x)dx 

l^ifci 0*0 1 l^ifca (x)\dx. 



Now we can use a classical argument ([35j,[36j[4]) to show that, for any r > 2, there exists 
C T > such that: 



L 2 (§ 2 ) 



\4>jkx ( x ) H^jfal i x ) dx 



< K R2j f 1 * 7 

" T U (1 + Bid (x, Q kl )) T (1 + BH (x, i jk2 )) T 



In order to evaluate this integral, we can for instance follow Q35J), by splitting the sphere 
S 2 into two regions: 

51 = {x e S 2 : d(x,£ jkl ) > d(fjfci,fjjfe 2 ) /2} 

5 2 = {x e S 2 : d(x,€ jk2 ) >d(£ jkl ,£ jha ) 1 2} . 
For what concerns the integral on Si, we obtain: 



Si 



1 + Bid (x, £ jkl )) T (1 + Bid (x, i jk2 )) T " (1 + Bid (£ ifclJ & fe )) T J Sl (1 + Bid (x, £ jk2 )) 



One also has that 



dx 



< 



Si 



1 + B>d (x, £ jk2 )) T ~ J § 2 (1 + Bid (x, £ jk2 )) 



dx 



2tt 



sin?? 



1 + Bi&) 



< 



< 



2vr 

2nC 
~&i 



V dy<** 



[i + y) 



£ 2 i 



ydy+ y T dy 



< 
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Because calculations on the region S 2 are exactly the same and because § 2 C Si U S 2 , we 
have that, for some constant C T depending on r, 



C 



+ BJd(t jkl ,£ jk2 )y 

yielding the desired conclusion. ■ 



Remark 5.2 Assuming that d (£jfc i; <v/fc 2 ) > $ uniformly for all j, we have immediately 

\ E $*X\<<,t 2 ><B- jT , 
where the constant k' t ^ 2 only depend on r, ( 2 ■ 

Remark 5.3 The previous Lemma provides a tight bound, of some independent interest, 
on the high frequency behaviour of covariances among wavelet coefficients for Poisson ran- 
dom fields. For Gaussian isotropic random fields, analogous results were provided by f&j/, 
in the case of standard needlets (bounded support), and by f2$f -F W{ . in the "Mexican" case 
where support may be unbounded in multipole space. It should be noted how asymptotic un- 
correlation holds in much greater generality for Poisson random fields than for Gaussian 
field: indeed in the latter case a regular variation condition had to be imposed on the tail 
behaviour of the angular power spectrum, and in the Mexican case this condition had to be 
strengthened imposing an upper bound on the decay of the spectrum itself. The reason for 
such discrepancy is easily understood: for Poisson random fields, non overlapping regions 
are independent, whence (heuristically) localization in pixel space is sufficient to ensure 
asymptotic uncorrelation; on the contrary, in the Gaussian isotropic case different regions 
of the field are correlated at any angular distance, and asymptotic uncorrelation for the 
coefficients requires a much more delicate cancellation argument. 



5.2 Fixed dimension 

Fix d > 2 and j > 1, consider a fixed number of sampling points {C^, £jfc d }, and define 
the associated (i-dimensional vector 

Pj- ■- {Pjki ' -iPjk d j> 

whose covariance matrix will be denoted by T t (note that, by construction, T t (i, i) = 1 for 
every i = l,...,d). Our aim is to apply the rough bound (12.1 7p in order to estimate the 
distance between the law of /?•. and the law of a random Gaussian vector Z ~ A/rf(0, If), 
where C = Id stands for the identity dx d matrix. Using Lemma IBTT] one has the following 
basic estimates: 
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c 



-II 



I op 



\c\\l = i 



\ fc 1 ^fc 2 =i 



— rf SU P 7i D /, ? — \yi 

fc 1 ^fc 2 =l,...,d Cjk 1 C r jk2 (.1 + -D- 7 " \Sjk 1 ,tLjk 2 )) 



< 



d 



X 



ac 2 



Ci?2 C 1 + 5i inffe 1 ^ 2 =i ) ... ! dC?(6fen 6*2)) 



Applying (12 . 1 T[) yields therefore that 

v ' 8 t^i ^ 2 



Mf (^) f{x)dx 



< A{t) + 



^I'M 8 



^J'fc llx, 3 (S 2 ) 



< A(t) + „,; 2 — BK 



where we used (13.211) and (I3.24p to yield <r| fc > Ci //2 ^ r l • We write this result as a separate 
statement. 



(5.29) 



Proposition 5.4 Under the above notation and assumptions, 



d 2 [pf t) ,Z\ < 



dC T ( 2 B^ T 



G<?2 (1 + i n ffc 1 ^fc 2 =l,...,d^(0fci'0 /c 2)) 



_ | (qi) 3 d% 2 y/2^ ? 



Because r can be chosen arbitrarily large, it is immediately seen that the leading term 
in the d 2 distance is decaying with the same rate as in the univariate case, e.g. B° j ' \[Rf 
Assuming however that d = d t , i.e. the case where the number of coefficients is itself 
growing with t, the previous bound may become too large to be applicable. We shall hence 
try to establish a tighter bound, as detailed in the next section. 
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5.3 Growing dimension 

In this section we allow for a growing number of coefficients to be evaluated simultaneously, 
and investigate the bounds that can be obtained under these circumstances. More precisely, 
we are now focussing on 

P j(t> \ p j(t)k^ "•>Pj(t)k dt )> 

where d t — > oo, as t — > oo. Throughout the sequel, we shall assume that the points at 
which these coefficients are evaluated satisfy the condition: 

„ , d te(*)fci^iW*a) ~ -JT ■ ( 5 - 30 ) 

Condition (I5.30p is rather minimal; in fact, the cubature points for a standard needlet /wavelet 
construction can be taken to form a maximal (c4)~ 1//2 -net (see |H [131 G51 SI] for more details 
and discussion). The following result is the main achievement of the paper. 

Theorem 5.5 Let the previous assumptions and notation prevail. Then for all t — 2,3..., 

there exist positive constants c and d , (depending on r, d, £ 2 but not from t,j(t), d(t) ) such 
that we have 

x(R t) „\ , cd t V2^ ddtB*® 



d 2 l3 y Z\ Z < , + — % - . (5.31) 

Proof. In view of (I2.16p and (15. 29ft . we just have to prove that the quantity 



8 ^ 3 / 2 «3 /p3 

is smaller than the second summand on the RHS of (15.311) . Now note that 

dt „ „ ( d t 



3 



J2 / l^Wi ( z )\ IV'iWfca ( z ) \ |^i(t)fe 3 (z)\dz<J2 \Y] (z)\ \ dz 

k lk2 kj^ 2 a JB(^ mx ,B-m) { k J 



where, for any z E B{^jU)\,B 



{i + Bmd{i mk ,z)} T 



< a^> + e 



.(0) {^' W %W*^(*)a)} T 
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Now for i- B(£j(t)\,B j( - t) ), x e B(£j(t)k,B we have by triangle inequality 

d (£j(t)k,£j(t)\) + d{£ j{t)k ,x) > d(£j(t)\,x), 
and because 

d (Z,j(t)k,£j(t)\) > d(£j( t )k,x), and 2d(Zj(t)h> £j(t)\) > d(£ j( t)\,x) , 
we obtain 

f ! / ^ ^ 

< V 1 f k t tbM dx<! J B m 

arguing as in [3], Lemma 6. Hence 

dt 

(*)| <<B i(t) , (5-32) 

it 

uniformly over z 6 S 2 , which immediately provides the bound. 

E / {x>iw*(*)l} ^<«^') 3 E / cfe = (^w) a . 

A JB{t mx ,B-j) { k J A JB(t mx ,B-m) 

Finally, to establish the sharper constraint 

J^fel^wOOlJ dz<K T d t B^\ 

it is sufficient to note that, exploiting ( 15.32p 

E /, ( 2 )| E l^'W*a (*)| E l^'(*)*s d * 

fci " /s2 fc 2 fc 3 

< ^B^Y^f \iP mkl (z)\dz = K 2 BWd m \\iP 
IT Js 2 



fcl 

< d t K 2 B 2 i®B-i® = d t K 2 B^ 



where we have used again H^jC^llLPfs 2 ) = 0(B 2 -'^ 2 p^ p ) = 0{B^ t ^ 2 - ) ), for p = 1. Thus 
(I5.3ip is established. ■ 

For definiteness, we shall also impose tighter conditions on the rate of growth of dt, B^ f ' 
with respect to R t , so that we can obtain a much more explicit bound, as follows: 
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Corollary 5.6 Let the previous assumptions and notation prevail, and assume moreover 
that there exists a, (3 such that, as t — >■ oo 

B 2j(t) ~R* ! <a< i ! d t &HP,0</3<a. 
There exists a constant k (depending on (1,(2, but not on j,dj,B) such that 

4 (iS5*),z)<«^, (5.33) 

for all vectors , fij^f), such that 115. 30\) holds. 
Proof. It suffices to note that 

dt ^ = 0(B- Tm d] +T/2 ) 

djBm / Rtd r \ v 2N 



= o 

and 



R t d t _ R t +I3T _ r -«+tW-«)+1 _ „m fnr t" ^ 1 ~ a 



Remark 5.7 From ( f5. it follows that for Rt — 10 12 we can establish asymptotic joint 
Gaussianity for all sequences of coefficients {^jm^uy ••' ^j(t)fc d (t)) °f dimensions such that 

0(1) , 



e.g. we can even at multipoles in the order 

of = O(10 3 ) we might take around 10 s coefficients with the multivariate Gaussian 
approximation still holding. These arrays would not be sufficient for the map reconstruction 
at this scale, but would indeed provide a basis for joint multiple testing procedures as those 
described earlier. 



Remark 5.8 Assume that dt scales as B 2 ^^; loosely speaking, this corresponds to the sit- 
uation when one focusses on the whole set of coefficients corresponding to scale j, so that 
exact reconstruction for bandlimited functions with I = O(B^) is feasible. Under this re- 
quirement, however, the "covariance" term A(t), i.e. the first element on the right-hand 
side of i5.31\) . is no longer asymptotically negligible and the approximation with Gaus- 
sian independent variables cannot be expected to hold. The approximation may however 
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be implemented in terms of a Gaussian vector with dependent components. For the sec- 
ond term, convergence to zero when dju) ~ B 2 ^^ requires B 3 ^^ = o[y/R~t). In terms of 
astrophysical applications, for R t ~ 10 12 this implies that one can focus on scales until 
180°/Bi ~ 180°/10 2 ~ 2°; this is close to the resolution level considered for ground-based 
Cosmic Rays experiments such as ARGO-YBJ (FTB{ ). Of course, this value is much lower 
than the factor B^ = o(^fR~ t ) = o(10 6 ) required for the Gaussian approximation to hold 
in the one- dimensional case (e.g., on a univariate sequence of coefficients, for instance 
corresponding to a single location on the sphere). 

Remark 5.9 As mentioned in the introduction, in this paper we decided to focus on a spe- 
cific framework (spherical Poisson fields), which we believe of interest from the theoretical 
and the applied point of view. It is readily verified, however, how our results continue to 
hold with trivial modifications in a much greater span of circumstances, indeed in some 
cases with simpler proofs. Assume for instance we observe a sample of i.i.d. random vari- 
ables {X t } , with probability density function /(.) which is bounded and has support in 
[a, b] C K. Consider the kernel estimates 

1 n Y — 

where K(.) denotes a compactly supported and bounded kernel satisfying standard regularity 
conditions, and for each j the evaluation points (x n o, x n Bj) form a B~~ J -net; for instance 

b — a 

a = x n0 < x nl ... < x nBJ = b , x nk = a + k ^~, k = 0, 1, B J . 



As argued earlier, conditionally on N t ([a,b]) = n, \5.3$ has the same distribution as 



L f Tyr I U Xnk > 



where N t is a Poisson measure governed by R t x J A f(x)dx for all A C [a, b\. Considering 
that |r — > a .s. 1) a bound analogous to 15. 331 can be established with little efforts for the vec- 
tor f n {x n ) '■= jin(^ni), •••> fn{ x nBi)j ■ We leave this and related developments for further 
research. 
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